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Abstract 



We establish a blow-up criterion in terms of the upper bound of the density 
Pu( ' and temperature for the strong solution to 2D compressible viscous heat-conductive 

■<^ . flo'vvs. The initial vacuum is allo'wed. 

-)— » 

ci '. 1 Introduction 

This paper is concerned with a blow-up criterion for the two dimensional compressible 
^ . viscous heat-conductive flows in (0, T) x Q 

m 

)0. ( pt + dw{pu) = 0, 



O ; ( c^{{pe)t + dw{p9u)) - kA9 + P{p, e)dwu = f I Vm + Vm^P + A(divM 



{pu)t + div{pu (g) n) + VP{p, 9) = fiAu + (A + p)Vdivu, (1.1' 



together with the initial-boundary conditions 

>< ; {p,u,e) = {po,uo,eo) in Q, (1.2) 

c^ ; (u,9) = (0,0) on {0,T) x dn. (1.3) 

Here fi is a bounded smooth domain in M?, and {p,u,6) are the density, velocity and 
temperature of the fluid, respectively. For a perfect gas, the pressure is given by 

p{p,e) = Rpe, (1.4) 

where i? > is a generic gas constant. The viscous coefficients p and A are constants 
satisfying 

p>0, X + p>0. (1.5) 
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Finally, c^, > and k > are the specific heat at a constant volume and thermal conduc- 
tivity coefficient, respectively. 

Some of the previous works in this direction can be summarized as follows. In the 
absence of vacuum, Matsumura and Nishida [25j proved the global well-posedness of the 
classical solution to (ll.ip - fll.5p with the initial data close to an equilibrium state and 
Danchin [6l [7] considered similar problems in the framework of critical spaces. With the 
aid of the effective viscous flux F, Hoff |13j proved the global existence of weak solutions 
with less restrictions on the initial data. If vacuum is taken into account, the problem 
becomes more complicated. Recently, Wen and Zhu [30] obtained a unique global classical 
solution to the ID model with large initial data and vacuum. For spherically symmetric 
flow in the exterior domain without the origin, the global existence of a strong solution 
was obtained by Jiang [12]. To our best knowledge, the ffist attempt towards the existence 
of weak solutions to the full compressible Navier-Stokes equations in dimension N > 2 
is given by Feireisl [11] where he proved the global existence of the so-called variational 
solutions in the case of real gases. In addition, it is worth mentioning that, by using a new 
mathematical entropy equality, Bresch and Desjardins [2] got the global weak solutions 
to the Navier-Stokes equations for heat conducting fluids with density and temperature 
dependent viscosity. Certainly, for the isentropic case, the results are more satisfactory, 
see [ini [H [HI ED 1221 [23 and references therein. 

On the other hand, when the initial density is compactly supported, Xin [31] proved 
that a smooth solution will blow-up in finite time in the whole space, see also [26] for 
a more general blow-up result. Thus, it is interesting to investigate the mechanism of 
blow-up and the structure of possible singularities. Some progress for the isentropic flow 
can be found in [T^l [ISl HSl [23 [2B] and references therein. For the non-isentropic case. 
Fan and Jiang [8] proved the following blow-up criteria for the local strong solutions to 
( ll.ip -( 1T3|l in the case of two dimensions 



or 



lim sup {IIpIIloo, Up II^oo, 116*11^:00} + / (||p||,4/i,go + ||Vp||j:2 + ||u||£,,L)(it = 00, 

T^T* \0<t<T Jo 

lim ( sup {\\p\\l^,\\p''^\\l^,\\6\\loc} + (||p||vi/i,?o + ||Vp||^2)c/t ) = cx) 

T^T* \0<t<T Jo ) 



provided 2p > A, where T* < 00 is the maximal time of existence of a strong solution, 
go > 3 is a certain number, 3 < r < cxd with 2/s -|- 3/r = 1, and L'''°° = U'°^{yC) is 
the Lorentz space. If the domain is a periodic or unit square domain in M^, the blow-up 
criterion is refined by Jiang and Ou [2D| to be 



T->-T* 



lim / ||Vu||Loo(it = 00, (1.6) 



which coincides with the Beale-Kato-Majda criterion ^ for ideal incompressible flows. 
For the case that fi is a bounded domain in R^, Fan, Jiang and Ou |9j established a 
blow-up criterion with additional upper bound of the temperature 6 



lim I sup ||6^||l°° + / llVMlUoofit = 00, (1.7) 

T^T* \o<t<T ' ' 



provided 

7fi > A. (1.8) 

Very recently, under the same condition (11.81) . Sun, Wang and Zhang [29] obtained a 
blow-up criterion in terms of the upper bound of (p, p^^, 6) 



lim sup {||p|U°° + Up ||l°° + ||^IU°°} = oo. (1.9) 

T^T* \o<t<T J 

However, from the physical point of view, it is natural to expect that the solution 
does not blow-up under conditions on the non-appearance of the concentration of the 
temperature and the density. For this reason, the aim of the current paper is to remove 
the lower bound of the density p in (II. 9p . As pointed out in [29], without the lower bound 
of the density p, it is difficult to deal with the highly nonlinear terms |Vm + Vm"""!^, (divw)^ 
in the temperature equation. To overcome this difficulty, we put together the estimates 
of supo<t<T llVwIl^a, supo<i<r \\^/pu\\l2 + /o^ IIVmII^j and supo<i<T ||V0|||2 + j^ \\pO\\h 
(see Lemma [3I2II3. 41 below) and make full use the good term J^ ||p^|li2- This is the main 
ingredient of our proof. Besides, let us emphasize that, instead of 6*, we use Ot as the test 
function in the proof of Lemma 13.41 which enable us to rewrite 



/ [^|VM + VM^P + A(divM) 

Jo 2 



as 

d /"* p 



dt 



I [ttIVm + Vu^\^ + A(divM)^]6' + other terms, 
Jo 2 

and thus we can avoid using the lower bound of the density. Finally, we remark that in 
the process of our proof, we have to deal with the troublesome term J^ || Vm||^4. To this 
end, we use the decomposition of the velocity u = v + w introduced by Sun, Wang and 
Zhang in [27]. More precisely, let v solves the elliptic system 

pAt; + (A + p) Vdivi; = VP{p, 6) in Q, 
f = on dVl. 

then from the momentum equation ( 11.1^ 9 and (Jl.lOp . it is easy to see that w solves 

pAw + (A + p) Vdivw = p-it in 17, 
w = on dQ. 

Here w is an important quantity, whose divergence could be regarded as the effective 
viscous flux and thus possesses more regularity information than u. Based on this decom- 
position, our approach next is to use the Gagliardo-Nirenberg inequality and then close 
the estimates with the help of Gronwall's inequality. Unfortunately, if the dimension 
A^ = 3, this approach does not work anymore, and that is why we only consider the two 
dimensional case. 



We would like to give some notations which will be used throughout the paper. 
Notations: 

1. Throughout this paper, we denote by fi a bounded smooth domain in R^ and U = 

Hi = {v ^ H^iyt) : \\v\\h^ < oo and v\qq, = in the sence of trace}, by virtue of 
Poincare's inequality, we redefine \\v\\hi = ||Vf ||l2. 

'^- f = ft + u ■ V/ is the material derivative of /, and / = (/)' denotes the twice 
material derivative of /. 

3. / / = / fdx and / f — / fdxdt. 

J Jn Jo J Jo Jn 

Before proceeding any further, we recall that Cho and Kim [3] obtained the local 

strong solution to fll.ll) - (ll.5p with initial vacuum for the spacial dimension A^ = 3 (and Q 

need not be bounded). Now if fi is a bounded smooth domain in M^ with (m, 6) satisfying 

the Dirichlet boundary condition (11. 3p . it is not difficult to verify that Cho and Kim's 

proof in [S] still works, and we state the corresponding result below. 

Proposition 1.1. Let q G (2, oo) be a fixed constant. Assume that the initial data satisfy 

Po>o, poeW^''^, iuo,9o)eHlnH^, 

and the compatibility conditions 

1 
/iAwo + (A + ii)Vdivuo - RV{po9o) = PoQi, 

kAOo + f |Vuo + VU([|2 + X{d%moy = p|^2, (1-12) 

for some {gi, g2) G L^ ■ Then there exist a positive constant Tq and a unique strong solution 
{p,u,9) to U.1]) - /[T75\) such that 

P>0, pGC([0,To];iyi'^), 
{u,e)eC{[0,To];HlnH^)nL\0,To;W^'i), 
{ut, Ot) G L2([0, To]; Hi), {^Ut, ^9^) G L-([0, To]; L^). (1.13) 

Our main result is stated as follows. 

Theorem 1.2 (Blow-up Criterion). Suppose that the assumptions in Proposition \1.1\ 
are satisfied, and assume that (p, u, 9) is the strong solution constructed in Proposition 
\l.li Let T* be the maximal existence time. If T* < oo, then 

lim sup (||p||l°o + ||^||l°°) = oo. (1-14) 



T^T 



0<t<T 



Remark 1.3. If either the boundary condition (11.31) on 9 is replaced by the Neumann 
boundary condition |^|af7 = or the bounded domain Vt is replaced by the whole space 
M^, we are able to obtain the same blow-up criterion as (I1.14p provided the corresponding 
local strong solution exists as in Proposition ll.il In these cases, we can not use Poincare's 
inequality directly to get the lower order estimates of u and 9, however, these would be 
recovered by using Lions's (for Vt bounded) and Hoff's (for Vt = M^) technics, and the 
readers are referred to [211 Remark 5.1 and fTH Lemma 2.3, for instance. 



2 Preliminaries 

Obviously, the appearance of the systems f ll.lOp and (11. lip make the estimates on the 
strongly elliptic operator fiA + (A + /i) Vdiv necessary. Now let U solves 

' fiAU + {X + n)VdiYU = F in Q, 
t/ = on 91]. 

First of all, we state some classical estimates for the above strongly elliptic systems, which 
will be used later frequently . 

Lemma 2.1 ([27J). Let p G (l,oo) and U be a solution of \2. 1\) . then there exists a 
constant C depending only on X,^,p and Q such that the following estimates hold: 

(1) ifFeLP, then 

\\U\\w'2,p<C\\F\\lv, (2.2) 

(2) ifF = divf with f = {fij)2x2, fij e LP, then 

\\U\\w^^r><C\\f\\LP. (2.3) 

As in [271 ESI El], we also need an endpoint estimate for the strongly elliptic operator 
/iA + (A+/i)Vdiv for the case p = oo. This will be done with the aid of the John-Nirenberg 
space of bounded mean oscillation whose norm is defined by 

ll/llBA^-OCf^) '■= ||/||l2{Q) + [f]BMO(n), 

with 

[f]BMO(n)-= sup I . ., / \f{y)- fnr{x)\dy, 
xen, re{o,d) l"ria^J| Jurix) 

fnrix) = TPTT^ / f{y)dy, 
I"r-I2;j| Jnr(x) 

where Qr{x) = Br{x) fl fi, Br{x) is a ball with center x and radius r, d is the diameter of 
Q and |r2j,(a;)| denotes the Lebesgue measure of Qr{x)- We have 

Lemma 2.2 (^). If F = divf with f = {fij)2x2, fij e L'^nL'^, then VU G BMO{n) 
and there exists a constant C depending only on A, /i such that 

\\VU\\BMOin)<C{\\f\\Loo + \\f\\L2). (2.4) 

Next, we state a variant of the Brezis-Waigner inequality [^ which, together with 
Lemma 12. 2[ will be used to give the gradient estimate of p. 

Lemma 2.3 ([27]). Let Q be a bounded Lipschitz domain in M^ and f G W^''^ with 
q G (2, oo), there exists a constant C depending on q such that 

II/IIl^ < C(l + ||/||BMoln(e + ll/lki,.)). (2.5) 

Remark 2.4. Lemma [2.21 and [2.31 are the two dimensional version of the corresponding 
lemmas in ^28j . To our knowledge, when the domain Q is bounded, the estimate (12. 4p can 
be found in [T] for a more general setting, and for the case fi = R^, Sun, Wang and Zhang 
gave a proof in [28] . The Brezis-Waigner type inequality (l2.5p was first established in 
on M^ and the case f2 be a bounded domain in M^(A^ = 2, 3) can be found in [271 [28] . 



3 Regularity of the velocity and temperature 

Let < T < T* be arbitrary but fixed. In what follows, we assume that (p, u, 6) is a 
strong solution of f ll.ip - fll.5p on [0,T] x Q with the regularity stated in Proposition 11.11 
Suppose that T* < oo. We will prove Theorem 11.21 by a contradiction argument. To this 
end, we suppose that for any T < T*, 

sup (IIpIU- + II^IU-) < M, (3.1) 

o<t<r 

where M is independent of T. We will deduce a contradiction to the maximality of T*. 
Hereinafter, we denote by C a general positive constant which may depend on the initial 
data, the domain Q, M in f l3.ip and the maximal existence time T*. 



rp 

sup [{p\u\^ + p\9\^)+ [ f{\Vu\^ + \V9\^)<C. (3.2) 

3<t<T J Jo J 



Lemma 3.1. Under assumption i\3. 1\) . there holds for any T < T* 

rp 

I {p\u\' + p\9\') + I 

0<t<T . 

The proof is the same as Lemma 2 in [29j , we omit it here. 
Lemma 3.2. Under assumption i\3. 1\) . there holds for any t < T* 

[\Vu\^+ [ fp\u\^<C + C f I ^\e\\Vu\+C f f \Vu\\ (3.3) 

J Jo J Jo J Jo J 

Proof. First of all, the momentum equation can be rewritten as 

pii + VP = pAu + (A + /i) Vdivu. (3.4) 

Taking the L^ inner product of the above equation with u, and integrating by parts, we 
get 

/ p|m|2 = / PdivM + / [pAu + (A + p)Vdivu\u. (3.5) 

Direct calculation yields 

/ PdivM = -f I i^divM - I {Pt + div(PM))divM 

+ I PVu: Vm^, (3.6) 

and 



[pAu + (A + p)Vdivn]'u 
--— / [p\Vu\^ + (A + p){dwuf] - p Vu: (VuVu) 
+ 1 f diYu\Vu\^ - {X + p) f diYuVu : Vu^ + ^^^ f{diYu)\ (3.7) 



Substituting fl3.6p and (13 .Tp into (13. 5p . we have 

j P\u\' + ~ y"[/i| Vm|2 + (A + /i)(divM)2] 
= — / PdivM - iPt + div(PM))divM + / PVu : Vu^ 
-jj Vu: (VnVn) + - / divM|Vup 
-(A + fi) divnVn : Vu^ + ^-^ / {divuf. (3i 



Notice that 

Pt + dw{Pu) = Rpd, (3.9) 

in view of the state equation P{p, 6) = Rp9 and (ll.ip i. Substituting (13. 9p into (13. Sp . and 
integrating the resulting equation over [0,t], we get 

/ /pN'+^/[Hv«r+(A+/^)(div«f] 

< C+ fpdivu + cf f y/p\9\\Vu\ + C f f\Vu\^ + C f f \Vu\^ 

J Jo J Jo J Jo J 

< C, + e [\Vu\^ + C [ I ^\0\\Vu\ + C f I \Vu\\ (3.10) 



where we have used Cauchy's inequahty, (13. ip and (13. 2p . After the second term on the 
right hand side of (I3.10p absorbed by the left hand side for a fixed e small enough, we 
complete the proof of Lemma 13.21 D 

Lemma 3.3. Under assumption Ii3. 1\) . there holds for any t < T* 

p\u\^+ [ [ \Vu\^ <C + C [ [p\9\^ + c[ [ \Vu\\ (3.11) 



Proof. Here the calculations are motivated by Hoff |J^. Taking the material derivative 
to (13. 4p . one deduces that 

pil + VPt + div(VP (g) u) 

= p[Aut + div(AM (g) u)] + (A + /i)[VdivMt + div((VdivM) m)], (3.12) 

where div(/ ® m) := ^ dj{fu^) as in |28l [29]. Taking the L^ inner product of the above 
equation with it, and using (ll.ip i. we arrive at 

-— / p|Mp = / Ptdivn- / div(VP®M)u 

+ / p[Aut + div(AM (g) u)]u +{\ + p) / [VdivMi + div((VdivM) (g) u)]u, (3.13) 



A short computation shows that 

/ Pidivn - / div(VP ^ u)u = / [P* + div(Pn)]divu - f PVu : Vm^, (3.14) 



nl^Ut + div(AM (g) u)\u 

= -/i / |VM|^ + /i /[(VmVm) : Vu + Vu : (VmVm) -divwVn : Vn], (3.15) 
and 

(A + /i) / [VdivMf + div((VdivM) (g) m)]m 

= -(A + /i) / IdivMp + (A + /i) / divwVM : Vii^ 

+ (A + /i) / [Vn : Vw^divu - (divu)MivM], (3.16) 

Combining (I3.13p - (l3.16p . and using (13. 9 p and (13. ip once more, we are led to 
--T- / p|m|^ + /i / |Vu|^ + (A + /i) / |divu|2 



= R p9divu - / PVm : Vm ' 

+/i / [(VmVu) : Vm + Vm : (VnVn) - divnVw : Vm] 

+ (A + yu) / [divwVM : Vm^ + Vu : Vu^divu — (divu)^div'u] 

< e / iVwp + a /" v^^' + a / iVwp + C, / |Vm|^ (3.17) 



Choosing a e sufficiently small, the ffist term on the right hand side of (I3.17P can be 
absorbed by the the left hand side, and then integrating the resulting equation over [0, t], 
using (13. 2 p , we complete the proof of Lemma 13.31 D 

Lemma 3.4 (Crucial estimates). Under assumption 1^3. 1\) . there holds for any T < T* 
sup [{\Vu\^ + p\u\^+\V9\^)+ [ l\Vu\^+l fp9^<C. (3.18) 

0<t<T J Jo J Jo J 

Proof. According to (ll.ip i. the thermal energy equation (II. ip ^:; can be rewritten as 

c^pO - kAO + RpOdiYu = ^\Vu + Vm^P + A(divM)2. (3.19) 

Multiplying the above equation by 9t, and integrating the resulting equation over Q, we 
obtain 



c^ pOu-Ve- R I p9dwu9 + R j p9divuu ■ VO 

^ /r/^IVT,. I VT, T|2 I x/i- „ \2i 



- / [/i(VM + Vm^) : (Vwt + VuJ) + 2AdivMdivMt]^. (3.20) 

Note that 

— 2A / divudivutO 

= — 2A / divudivuO + 2A / divMdiv(M ■ Vu)9 

= — 2A / divudivuO + 2A / divudiU^dku'-Q + 2A / divu-u''(9fcdivu^ 

= -2A / divwdivM^ + 2A / divudiu'^dku'e + A / ^m ■ V(divM)2 

= -2A / divudiviiO + 2A / divuVu : Vm^^ - A / ^(divw)^ - \ u- V^(divM)^ 

(3.21) 

and 

- /i / (Vm + Vm"^) : (Vwi + VuJ)e = -2/i / Vm : VMi6' - 2/i / Vm : Vm7^ 

= /1 + /2, (3.22) 



Ji = -2/i / Vm : Vm^ + 2/i / diuWi{u''dkU^)e 

= -2/i / Vm : Vm6' + 2p Vu: {VuVu)9 + p j Ou ■ V| Vm|^ 

= -2/i / Vm : Vm^ + 2/i / Vm : (VmVm)6' - /i / MivM] Vwp - p u- V6'| Vm^, 

I2 = -2/i /" Vm : Vii'^e + 2p f diU^dj{u^dku')e 

= -2/i f Vu : Vm^^ + 2/i /" diU^dju'^dku'e 

= -2/i / Vm : Vm^^ + 2/i [vu: {VuVufO + p [ Ou^dkidiU^dju') 



-/i / edivuVu : Vu^ - f^ u ■ VOVu : Vn^. (3.23) 



Substituting fl3.2ip - fl3.23p into f l3.20p . using Cauchy's inequality and f l3.2p . we get 

Cv l _h2 , n I IV7„,|2 , r^ / L,|2|v7fl|2 , ^ / |V7„M2 



< -^ I p9' + C, I \Vu\' + C I |M|^|V^r + e / |Vm 



+C I |Vm|^ + J [[^\Vu + Vm^P + X{diYuf]e. (3.24) 



Consequently, integrating f l3.24p over [0,t] and using (13. 2p again, we are led to 
Jo J J Jo J Jo J 



+C I f\Vu\^ + C j \Vu\\ 
Substituting (13. 3p and (13. lip into the above inequality, we have 

' I pO^+ /iv^p 

< C, + C I f \u\^\V9\^ + eC f f p\9\^ 

+c I I yfp\e\\Wu\ + c I I \wu\^ 
< c, + c [ [ \u\^\ve\^ + ec [ [p\e\^ + c [ I \vu\\ 

Jo J Jo J Jo J 

where we have used Cauchy's inequality and (13. 2p . Take a e sufficiently small, then there 
holds 

2 I / iv7Z)|2 ^ r^ t r^ I / i„.|2iv7Z)|2 , r^ I iv7„.|4 



pr+/|v^r < c + c \u\'\ve\' + c \Vu\\ (3.25) 

J J Jo J Jo J 

It follows from Lemma 12.11 Gagliardo-Nirenberg' inequality and (13. ip that 

< C 



[ f\P\^ + C [ ||Vw;||i2||Vw||^i 
Jo J Jo 

+ C [\\\Vu\\l. + \\Vv\\l.)\\^u\\l 
Jo 



10 



In particular, 



Moreover, 



< C + C f\\\P\\l, + \\Vu\\h)\\^pu\\l. 

Jo 

< C + C [ il + \\Vu\\l,)\\^u\\l,, (3.26) 

Jo 



\VU\\% < C\\P\\l, + Ci\\P\\L^ + \\Vu\\L^)\\^u\\L2 

< C + Cil + \\Vu\\L2)\\^u\\L^. (3.27) 



t p pt pt pt 

u\^\V9\^ < / ||n||ioo||V^||i2 <C / \\u\\l^+C ||V^||i2. (3.28) 
J Jo Jo Jo 



By Sobolev imbedding W^'^ ■— )■ L'^ 



\u\ 



< C\\u\\% + C\\Vu\\% < C\\Vu\\%, (3.29) 



where we have used Poincare's inequahty. Then (I3.28P and (I3.29P imply that 

f f\u\^\V9\^ < C [ {\\Vu\\l, + \\V9\\i,). (3.30) 

Jo J Jo 

Substituting (13.261) and (13.301) into (13.251) . we obtain 

< C + C [ \\ve\\l, + c [ (1 + liv^iiioil v^^'lli- (3.31) 

Jo Jo 

Now we are in a position to close all the estimates above. Indeed, 4Cx (l3.3ip + ( lX^ + (l3.1ip 
implies 

1 2 I „l • |2 I \T-7n\2\ I / / |V7-, |2 I / / „n2 



< C + C f \\V0\\l + C [\l + \\Vu\\l,)\\^u\\l, 

Jo Jo 

< C + C f (1 + \\Vu\\l, + \\V9\\h){\\Vu\\l, + \\,/pu\\h + llV^lliO. (3.32) 

Jo 

By virtue of Gronwall's inequality and (13. 2p . we complete the proof of Lemma [3.41 D 
Corollary 3.5. Under assumption jiS. 1\) . there holds for any T < T* 

sup (||m||l°c + ||Vm||l4) < C, (3.33) 

0<t<T 

I \V^9\^ < C (3.34) 

J 

11 



Proof. Clearly, fl3.33p is a direct consequence of fIS.lSp . fl3.27p and fl3.29p . In view of 
( I3.18p . fl3.33p and f l3.ip . applying the standard elliptic regularity theory to (13.1 9p . then 
( I3.34P follows immediately. D 

Lemma 3.6. Under assumption liS. 1\) . there holds for any T < T* 

sup [pO^+f /|V^P<C. (3.35) 

0<t<T J Jo J 

sup ||^||j:^2 < C. (3.36) 

0<t<T 

Proof. Taking the material derivative to (I3.19p . we find 

c^p9 - k{A9)- = [cy - R)pmNu + Rpe\{^\^uf + Vm : Vm"^] 

- RpMwii + [/i(VM + Vm^) : ( Vm + Vm^) + 2AdivMdivM] 

- [/i(VM + Vm^) : (VmVm + Vm^Vm^) + 2AdivMVM : Vm^], (3.37) 
where we have used the facts 

p = —pdivu, (div-u)' = divit — Vm : Vu^ 

and 

(Vm)' = Vn - VwVn, {Vu^J = Vii^ - Vn^Vw^. 

Taking the L^ inner product of the (I3.37P with 6, and using (ll.ip i. we arrive at 

< (c^ -R) f pe^ivu + R I pee[{divuf + Vu ■. Vu^] 

-R f pOedivu + C f |Vm||Vm||^| + C f \Vu\^\e\. (3.38) 

-K HAeye = -k nAOt + u ■ vA9)e 

K I VOtVe + K [ diU^dkiOe + k [ u^dkiOdiO 



Note that 



K I \V9\^ -K I vevuve + k j vu- v^ee. (3.39) 



We infer from (EJHD, (EJSI) and (^ that 



fjj^' + '^im 



12 



< C f p9^\Vu\+C f p\9\\Vu\^ + C f p\9\\Vu\+C f \Vu\\Vu\\9\ 
+C I \Vu\^\9\ + K I |V^||Vm||V^| + k f |Vm||V^^||^| 

< f / I V^P + e f \Vu\^\9\^ + cJ p9' + cj I Vm|^ 

+C, f |Vm|2 + C, f \V^9\^ + C f \Vu\^\V9\\ (3.40) 

From (13.331) . Gagliardo-Nirenberg' inequality and Poincare's inequalty, we have 

e f \Vu\^\9\^ < e||Vn||i4||^||i4 <eC||^|U2||V^||L2 <C||V^||i2, (3.41) 

and 

/"iVmHV^I^ < \\Vu\\14V9\\1,<C\\V9\\l2\\V9\\hi 

< C||v^||i2 + c||v'^||i2<c + c||v'^lli2, (3.42) 

where we have used (I3.18P in the last inequality above. Substituting (13.411) and (13.421) 
into (I3.40p . choosing a e small enough and using (13.331) . we obtian 

c,^ f p9^ + K f \V9\^ <C + C f \V^9\' + C f iVul' + C f p9\ 

Integrating the above equation over [0,t], and using (I3.18P and (I3.34p once more, we get 
(13:351) . Finally, (1336|) follows from (jSHSD, (I3A8|) . ([233]), (l3:35|) and ^^ immediately. 
This completes the proof of Lemma 13. 6[ D 

4 Improved regularity of the density and velocity 

Generally speaking, for the continuity equation (ll.ip i. the gradient estimate on p rely 
on the boundedness of /^ ||Vm||lc=o. On the other hand, due to the lower regularity on 
p, we can not obtain the higher order regularity of u through the momentum equation 
(13. 4p . Nevertheless, with the help of the decomposition m = t> + w, we can close the 
gradient estimate of p based on a logarithmic estimate for the strongly elliptic operator 
/iA + (A + /i) Vdiv. 

We first give some higher regularity on w below. 



sup||w||h2+/ (||V^w||i, + llVwIlioo) < C, gG(2,oo). (4.1^ 

3<t<T Jo 



Lemma 4.1. Under assumption \3. 1\) . there holds for any T < T* 

0<t<T Jo 

Proof. Using (12. 2p with U = w,F = pii and p = 2 , we have 

\\w\\h'2<C\\pu\\l2<C\\^u\\l^, (4.2) 
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where we have used (13.11) . Then it follows from fl3.18p that 

sup ||w||h2 < C. (4.3) 

0<t<T 

For q > 2, using (12. 2p . (13. ip again, and by virtue of Gagliardo-Nirenberg inequality, we 
find 

||V^w||l? < C||pm||l9 < C||m||l, 

< C\\u\\%\\Vu\\7V < C\\Vu\\l2, 

where we have used Poincare's inequality. Next, we use the Sobolev imbeding W^''^ M- L°° 
for g > 2 to get 

\\Vw\\l^ < C{\\Vw\\l. + WVMlo) < C{\\Vw\\h^ + \\VMl'>), (4.4) 

then we infer from fl318|) . K3^ and (531) that 



/ i\\^M\h + \\^w\\U)<c. 

Jo 



(4.5) 



This completes the proof of Lemma 14. 1[ D 

Thanks to all the estimates obtained above, we will get the gradient estimates of the 
density p next. 

Lemma 4.2. Under assumption ^3.1\) . there holds for any T <T* 

sup (||p|ki>« + \\u\\h^) + / llVwIlioo < C. (4.6) 

0<t<T Jo 



Proof. The proof follows the ideas of Sun, Wang and Zhang [271 128| [29] , we sketch it here 
for completeness. First of all, take V to the continuity equation (ll.ip i to find 

dtVp + {u ■ V) Vp + VwVp + divwVp + pVdivw = 0. (4.7) 

Multiplying (14. 7p by g|Vp|''^^Vp and integrating over Q, we obtain 

^||Vp|U, < C\\Vu\\l^\\Vp\\l.+C\\\/\\\l. 

< C{\\Vv\\l^ + II V«;|Uo.)||Vp|U, + C{\\VMl^ + W'^Mlo) 

< C(l + ||Vt;||Loo + ||Vw||loo)||Vp||l, + C(||VV»||l, + ||V0||lO, 

(4.8) 

where we have used (12.20 . To close (14. 8p . we have to bound || Vf ||l°o, and it is just this term 
leads us to show the endpoint estimate for the strongly elliptic operator pA + (A + p)Vdiv. 
In fact. Lemma [2. 1112.31 imply that if g > 2 

||Vt;|Uo. < C(l + ||Vt;||BA/oln(e+||V\||LO) 
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< C(l + (llPlUoo + \\P\\l2) ln(e + \\VMl^)) 

< C(l + ln(e + II VplUO + ln(e + || V^|U,)) 

< C(1 + ||V0|U. +ln(e+||Vp|UO)- (4-9) 
Substituting iKdh into (H^ . we get 

^(e + llVplUO < C(l + IIV^IU. + ||Vw|Uoo)||Vp|U, + Cln(e + ||Vp||l,)I|Vp||l. 

+C{\\V'w\\l'> + \\'^0\\l,), (4.10) 

then using (14.11) . (I3.36P and Gronwall's inequality, one deduces that 

sup ||Vp|U, < C. (4.11) 

0<t<T 

From (13:361) . fHlD . fHTTj) and Lemma EH there holds for g > 2 

iVwIlioo < C ||Vt;||ioo + / ||Vu;||ioo 
Jo Jo 

|2 , ry IIV72„.||2 



< c + c \\vv\\i, + c jiv'ujii, <a (4.12) 

Jo Jo 

Moreover, it follows from (O), (Km and fHTTl) 

||V^n||i2 < ||V^^||l2 + \\V\\\l^ <C + C||Vp|U2 < C. (4.13) 

This completes the proof of Lemma 14. 2[ D 



5 Proof of Theorem 11.2 



The combination of Lemma 13. 4[ Lemma 13.61 and Lemma 14.21 will enable us to extend the 
strong solution (p, u, 6) beyond the maximal existence time T*. Indeed, by virtue of (14. 6p . 
(13.361) and time continuity stated in (I1.13p . we can define 

(p,M,6')|t=T. = lim(p,u,6'), (5.1) 

and 

h := pu\t=T* = lim (pA-u + (A + p)VdivM — VP) strongly in L^. (5.2) 

On the other hand, using Sobolev imbedding W^''^ M- C ' ^9 for g > 2, (14. 6 p and the time 
continuity on p stated in (11.130 . one easily deduces that 

p|i=T*(a^) = lim p(x) uniformly in ^2, 



and hence 



^|j=2^*(x) = hm ^/p{x) pointwise in Q, (5.3) 
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then (15. 3p . the assumption (13. ip that p is upper bounded and Lebesgue's Dominated 
Convergence Theorem imply that 

^/p\t=T* = hm ^/p strongly in L?. (5.4) 

Besides, we infer from (I3.18P that there exists a sequence t^ — )■ T* as fc — )■ oo and a 
function gi E L^, such that 

gi = lim {^/pu){tk) weakly in L^. (5.5) 

It follows from flOl) and fl575l) that 

A/p|t=T*^i = lim {pu){tk) in the sence of distribution. (5.6) 

Comparing (15. 2 p with (15. 6p . we obtain 

h = ^/p\t=T*9u 

i.e. 

(/iAu + (A + /i) VdivM - VP) |i=T* = ^\t=T*~gi- (5.7) 

Similarly, there exists a function ^2 ^ -^^ such that 

{kM + ||Vm + Vm^|2 + A(divM)2)|,=3.. = y/-p\t=T'~g2- (5.8) 

Now (15. ip . (15. 7p and (15. 8p assure that we can take (p, m, 6')|j=t* as the initial data and 
apply Proposition 11.11 to extend the local strong solution beyond T*, which contradicts 
the maximality of T*. This completes the proof of Theorem 11.21 D 
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